HEWLETT-PACKARD 


HP-A1C 


MATHEMATIK-PAKET 


Hinweis 


Das hierin enthaltene Material ist mit keiner Verpflichtung oder Garantie 
irgendeiner Art verbunden. HEWLETT-PACKARD übernimmt infolge- 
dessen keine Verantwortung und wird keine daraus folgende oder 
sonstige Haftung übernehmen, die auf irgendeine Art aus der Benutzung 
dieser Programmsammlung oder Teilen davon entsteht. 


HEWLETT-PACKARD bittet um Ihre Aufmerksamkeit 


Um die Software-Unterstützung zu verbessern, bitten Sie die HP-Appli- 
kationsingenieure um Ihre Hilfe. Ihre frühzeitige Antwort hilft uns, die 
Qualität der Software und existierende Programm-Pakete zu verbessern. 
Ihre Beantwortung dieses Fragebogens ist aüßerst wertvoll für uns. 
1. Name des Paketes: Mathematik 
2. Wie bedeutsam war die Verfügbarkeit dieses Paketes für den Kauf 
eines Hewlett-Packard Rechners? 
U hätte ihn ohne dieses Programm-Paket nicht gekauft 
O wichtig DU unwichtig 
3. Was ist das Hauptanwendungsgebiet für dieses Paket? 


4. Geben Sie bitte in der untenstehenden Liste den von Ihnen empfun- 
denen Gebrauchswert der Programme an: 
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5. Haben Sie einen Drucker gekauft? Dja Dlnein 
Wenn ja, ist das Druckformat dieses Paketes zweckmäßigD jaU nein 
6. Welche Programme fehlen nach Ihrer Meinung in diesem Paket? 


7. Welche weiteren Programm-Pakete wünschen Sie sich? 


DANKE FÜR IHRE ZEIT UND MITARBEIT. 


Name Stellung 
Firma 
Straße (Land) 


PLZ, Stadt Telefon 


Zusätzliche Bemerkungen 


Hewlett-Packard SA 
PERSONAL CALCULATOR DIVISION 
P. O. Box 


CH 1217 Meyrin 2 -— Geneva 


EINLEITUNG 


Die Programme des Mathematik-Pakets stammen aus den Gebieten Funktio- 
nentheorie, numerische Verfahren, lineare Gleichungssysteme, analytische Geo- 
metrie und ausgewählte Funktionen. Jedes Programm dieses Pakets besteht aus 
einem eigenständigen Programm im Anwender-Modul und einem Abschnitt 
in dem vorliegenden Bedienungshandbuch. Das Bedienungshandbuch enthält 
die wichtigsten Gleichungen und die Anweisungen zur Benutzung des Pro- 
gramms. Außerdem werden eines oder mehrere durchgerechnete Beispiele mit 
der zur Lösung benötigten Liste von Tastenfolgen beigefügt. 


Schalten Sie bitte den Rechner dann ab, wenn Sie ein Anwender-Modul ein- 
stecken und vergewissern Sie sich, daß Sie mit dem Abschnitt „Einstecken 
und Entfernen der Anwender-Module“ vertraut sind. Wenn Sie ein Programm 
verwenden wollen, so sollten Sie sich die Zeitnehmen, die Abschnitte „Formatder 
Benutzer-Anweisungen“ und „Einige Anmerkungen zur Programm-Verwen- 
dung“ durchzulesen. 


Zunächst sollten Sie ein Programm dadurch kennenlernen, daß Sie es ein- 
oder zweimal anhand der vollständigen Liste von Benutzeranweisungen durch- 
laufen lassen. Danach genügen Ihnen sicher die Anzeigen des Programms 
oder die Kürzel auf den mitgelieferten Masken für das Tastenfeld, aus denen 
hervorgeht, welche Variablen einzugeben, welche Tasten zu drücken sind und 
welche Werte ausgegeben werden. Zur Erleichterung Ihrer Arbeit haben wir 
eine Karte mit Kurzbeschreibungen für jedes Programm - einschließlich der 
Benutzer-Anweisungen - beigelegt. 


Wir hoffen, daß Sie von dem Mathematik-Paket bei der Lösung zahlloser Pro- 
bleme ihres Fachgebietes unterstützt werden. Eine Bittenoch: Wirsind sehr daran 
interessiert, Ihre Bewertung der Programme dieses Pakets zu erfahren, weswegen 
wir einen Fragebogen vorbereitet haben, der sich in der vorderen Umschlag- 
decke des Bedienungshandbuches befindet. Dürfen wir Sie bitten, diesen 
Fragebogen mit Ihren Kommentaren zu diesen Programmen auszufüllen? Eben 
aus Ihren Kommentaren können wir lernen, wie wir die Brauchbarkeit und den 
Wert unserer Programme verbessern können. 


INHALT 
IE nen a ee nee ee ee 1 
a ee ERORFEREERFERRGRG SEE ERDE EEE TOTEN N TE ET 2 
Einstecken und Entfernen der Anwender-Module......eescsccserneenee 4 
Format der Benutzer-Anweisungen ..o.s-soososeeeeeeeerennnnnnn nenne 6 
Einige Anmerkungen zur Programm-Verwendung..csesrcseereonnnnnnee 8 
Malie Free 0 > 10 


Es werden Determinante und Inverse für Matrizen bis zu 14 Reihen und 
die Lösungen für simultane Gleichungen bis zu 14 Unbekannten be- 
stimmt. 


Linus ir d-Oiin ine ara... en anna 17 
Dieses Programm verwendeteinenabgekürzten NEWTON-Algorithmus, 
um die Lösungen einer gegebenen Gleichung zu finden. 


Nullstellen und Funktionswerte von Polynomen „os ceseseeeneennnonene 21 
Im Programm werden alle Nullstellen von Polynomen mit reellen K.oefli- 
zienten höchstens fünften Grades bestimmt. Es können auch Funk- 
tionswerte für beliebiges x berechnet werden. 


Ps a RE 5 ee nf ee we 25 
Das Programm leistet die numerische Integration sowohl für explizite 
Funktionen, als auch für eine endliche Anzahl von Punkten gleichen Ab- 
stands (diskrete Werte). Die Integrale expliziter Funktionen werden nach 
SIMPSON, diskrete Fälle entweder nach der Trapez-Regel oder nach der 
SIMPSONschen Regel bearbeitet. 


EEE nn an anna ae ee 29 
In diesem Programm werden Differentialgleichungen ersten und zweiten 
Grades nach der RUNGE-KUTTA-Methode gelöst. 


IIND: sen rer iin na ann 34 
Dieses Programm berechnet Fourier-Koeffizienten aus Stichproben 
periodischer Funktionen. Bis zu zehn Koeffizientenpaare können 
simultan aus N gleichabständigen Ordinatenwerten berechnet werden. 
Die Koeffizienten können sowohl als kartesische als auch als Polar- 
Koordinaten ausgegeben werden. 


BEE nn ri ee 38 
Diese Programmsammlung ist für konjugiert komplexe Zahlen bestimmt. 
Neben den vier Grundrechnungsarten (+,-,x,*) können die gebräuch- 
lichsten Funktionen der komplexen Variablen z und w bestimmt werden 
(wie: |z’], Vz, z", z”", e?, In z, sin z, cos z, tan z, a?, log.?, z”“ und z”). 


Inhalt 3 


isses: Tri == 2a ar re aan 44 
Dieses Programm berechnet die hyperbolischen Funktionen sinh x, 
cosh x, tanh x, sowie deren Inverse. 


iii = 5.55 5: 5 ee Tee ee 46 
Diese Programme bestimmen Fläche, Seitenlänge und Winkel für jedes 
definierte ebene Dreieck. 


se = 25 2 nn na en a a a ee 52 
Dieses Programm leistet 2- und 3- dimensionale Koordinaten-Translation 
und/oder Achsenrotation. 


WU TER ER ur nen 57 


EINSTECKEN UND ENTFERNEN DER ANWENDER-MODULE 


Machen Sie sich bitte mit der nachfolgenden Information vertraut, bevor Sie 
ein Anwender-Modul zum ersten Mal einstecken. 

Es können bis zu vier Anwender-Module in die Buchsen des HP-AIC eingesteckt 
werden. Sind diese Module eingesteckt, so können die Namen ihrer Programme 
durch das Drücken von U) 2 in die Anzeige abgerufen werden. 


VORSICHT 


Schalten Sie den HPA41C immer vor dem Einstecken oder Entfernen einer 
steckbaren Erweiterung (Anwender-Modul oder auch Zubehör) aus. Versäumen 
Sie das Ausschalten, so kann sowohl der Rechner als auch das Zubehör be- 
schädigt werden. 


Im folgenden ist beschrieben, wie Sie ihr An- 
wender-Modul einstecken müssen: 


1. Schalten Sie den HP-A1C aus! Versäumen 
Sie dies, so kann außer dem Modul auch 
der Rechner beschädigt werden. 


2. Entfernen Sie die Schutzkappen auf den 
Buchsen, Bewahren sie die Schutzkappen 
auf, da freie Buchsen durch diese Kappen 
verschlossen werden sollen. 


3. Stecken Sie das Anwender-Modul in eines 
der Buchsen, jedoch immer nach dem 
letzten verwendeten Speichererweiterungs- 
modul. 


Einstecken und Entfernen der Anwender-Module 


4. Wollen Sie weitere Anwender-Module ein- 
stecken, so können diese in jede beliebige 
Buchse eingesteckt werden, jedoch immer 
nach dem letzten Speichererweiterungs- 
Modul. Haben Sie beispielsweise ein Spei- 
chererweiterungs-Modul (RAM) in Buchse 
1 eingesteckt, so können Anwender- 
Module (ROM) nur in die Buchsen 2, 3 
und 4 beliebig eingegeben werden. . 
Stecken Sie niemals ein Anwender-Modulin 
eine Buchse mit niedrigerer Nummer, als 
ein Speichererweiterungsmodul. 
Vergewissern Sie sich, daß freie Buchsen 
durch Schutzkappen verschlossen sind. 


5. Schalten Sie den Rechnereinund folgen Sie 
den programmspezifischen Anweisungen 
dieses Handbuches. 


Entfernen der Anwender-Module: 

1. Schalten Sie den HP-41IC aus! Versäumen 
Sie dies, so können Rechner und Modul be- 
schädigt werden. 


2. Nehmen Sie das entsprechende Modul und 
ziehen Sie es heraus, wie auf der Abbil- 
dung dargestellt. 


3. Stecken Sie eine Schutzkappe auf die un- 
belegte Anschlußbuchse. 


MISCHEN VON SPEICHERERWEITERUNGS- (RAM) 


UND ANWENDER-MODULEN (ROM). 
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Die Auslegung des HP-41C erfordert es, daß RAM-Module in den Buchsen mit 
niedrigerer Nummer einzustecken sind. Dies gilt für alle Erweiterungen, auch für 


solche wie den HP-82104A Kartenleser oder den HP-82143A Drucker. 


Verwenden Sie Speichererweiterungs-Module und Anwender-Module, so müs- 
sen die Speichererweiterungs-Module immer in die Buchsen mit der niedrigeren 
Nummer eingeschoben werden, die Anwender-Module in die Buchsen nach dem 
letzten Speichererweiterungsmodul. Wollen Sie RAMs und ROMs gemischt ver- 
wenden, so können Buchsen freigelassen werden. Beispielsweise,stecktin Buchse 
lein RAM, so kann ein Anwender-Modul in Buchse 4 unter Auslassung der 


Buchsen 2 und 3 eingeschoben werden. 


FORMAT DER BENUTZER-ANWEISUNGEN 


Die vollständige Tabelle der Benutzer-Anweisungen - zu jedem Programm gibt 
es im Bedienungshandbuch eine Tabelle - enthält die Anleitung zur Bedienung 
eines jeden Programms aus diesem Paket. 


Die Tabelle besteht aus fünf beschrifteten Spalten. Von links nach rechts gelesen, 
enthält die erste Spalte die laufende Nummer der Anweisung. 


Die Spalte ANWEISUNG enthält die Anweisung und kommentiert die auszu- 
führenden Schritte, 


In der Spalte EINGABE sind die Eingabedaten spezifiziert, gegebenenfalls deren 
Maßeinheiten oder die entsprechende alphanumerische Reaktion auf eine Ihrer 
Antworten. Die Tasten für die Dateneingabe sind die Tasten für 0 bis 9 (nume- 
rische Tasten), [EEX]) (Eingabe Exponent) und (Vorzeichenwechse!). 


Die Spalte FUNKTION gibt die Taste(n) an, die nach der Eingabe der ent- 
sprechenden Daten zu drücken sind. 


Immer dann, wenn eine Angabe in den Spalten EINGABE oder FUNKTION 
goldfarben gedruckt ist, ist zuvor die [ALPHA ]-Taste zu drücken. Sobald die 
Informationen eingetastet sind, muß wiederum die [ALPHA J-Taste betätigt werden, 
um den Rechner in seinem normalen Modus umzuschalten, oder auch um die 
Programmausführung zu starten. Beispielsweise meint IUU N, daß die 
Tastenfolge [XEQJ[AtrHA JFOUR zu drücken ist. 


Die Spalte ANZEIGE gibt die Reaktionen des Rechners sowie Zwischen- und 
Endergebnisse mit ihren Maßeinheiten an, sofern dies zutrifft. 


Über der Spalte ANZEIGE steht ein Kasten, der die Mindestzahl von Registern 
(Speicherplätzen) enthält, die zur Programmausführung notwendig sind. Die 
Seiten 73 bis 117 des Bedienungs- und Programmierhandbuches enthalten 
detaillierte Informationen darüber, wie man die belegte Anzahl von Registern 
bestimmt. 


Im folgenden ist die Tabelle von Benutzer-Anweisungen des Programms 
FOURIER-Analyse dargestellt. 


Umfang: 027 


Programm initialisieren NO.SAMPLES=? 
Anzahl der Stichproben einer Folge 


eingeben NO.FREQ=? 
Anzahl gewünschter Frequenzen 

eingeben 1ST.COEFF=? 
Exponent des ersten Koeffizenten (J). Y1=? 

Eingabe von Y,„. n1, ... ‚N Y2=", ... ‚RECT? 
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Schritt 5 wiederholen bis RECT? in der 
Anzeige erscheint. 


Ist die Antwort ja, drücken der Taste 
um die Koeffizienten für 
JxSJ+ # fregs. in kartesischen 
Koordinaten auszugeben. 


ist die Antwort nein (N), werden die 
Koeffizienten in Polarform ausgegeben. 


a 
17) 


EI 
17210177] 


drücken von zeigt aufeinander- 
folgende Koeffizienten an. 


Um den Funktionswert für t zu berech- 
nen, USER-Modus setzen und t eingeben. 


u: 


EINIGE ANMERKUNGEN ZUR PROGRAMM-VERWENDUNG 


KATALOG 


Sobald ein Anwender-Modul in eine Buchse des HP41C eingesteckt ist, kann 
sein Inhalt durch Drücken von @EaAtauog 2 (der Katalog der Erweiterungen) 
ausgegeben werden. Drücken der Eatauogl-Funktion bewirkt, daß eine Liste 
aller Programme und Funktionen des Moduls aber auch aller Funktionen jeder 
anderen eingesteckten Erweiterung ausgegeben wird. 


MASKEN 


Für eine Anzahl der Programme des Pakets wurden Masken beigelegt. Um diese 
Programme zu benutzen, wählen Sie die entsprechende Maske und legen Sie sie 
auf den Rechner. Die Bezeichnungen auf den Masken sollen Ihnen bei der Ver- 
wendung des Programmes helfen. Der Programmname ist vertikal auf der linken 
Seite der Maske angegeben. Die Bezeichnungen in blau beziehen sich auf die 
unmittelbar über ihnen stehenden Tasten, wenn die Maske richtig plaziert ist 
und der Rechner im USER-Modus betrieben wird. Goldfarbene Bezeichnungen 
haben die gleiche Bedeutung, nur sind sie oberhalb der entsprechenden Taste 
angebracht, und zusätzlich ist vorher die Präfix-Taste zu drücken. Um es noch- 
mals zu betonen, der USER-Modus muß eingestellt sein. 


ALPHA- und USER-Modus 


In diesem Bedienungshandbuch ist eine besondere Notation verwendet, um den 
ALPHA-Modus zu kennzeichnen. Immer dann, wenn auf der Benutzer-An- 
weisungs-Tabelle eine Anweisung goldfarben gedruckt ist, muß vor der An- 
weisung die [ALPHA] -Taste gedrückt werden. Ist die Anweisung eingegeben, 
ist [AtPHA] nochmals zu drücken, um den Rechner in seinen normalen Be- 
triebszustand zurückzuversetzen oder eine Programmausführung zu starten. 
So bedeutet KEG) FOUR z. B., daß die Tasten FOUR zu 
drücken sind. 


Im USER-Modus, bezogen auf die neudefinierten Tasten und deren obere Be- 
schriftung, werden in diesem Bedienungshandbuch die Symbole [A] - und 
MlAl- WIE) in Benutzer-Anweisungs-Tabellen und bei gelösten numerischen 
Beispielen verwendet. 


VERWENDUNG DES ZUSÄTZLICHEN DRUCKERS 


Ist der als Zubehör erhältliche Drucker an den HP-AIC angeschlossen, so werden 
alle Ergebnisse des Mathematik-Pakets Anwender-Moduls automatisch aus- 
gedruckt. Wünschen Sie zusätzlich noch eine Liste aller Eingabewerte bei einem 
bestimmten Programm, so schieben Sie den Drucker-Wahlschalter in Stellung 
NORMAL, bevor Sie das Programm ausführen. In diesem Modus werden alle 
Eingaben und Tastenfolgen auf dem Drucker gelistet; somit verfügen Sie über 
eine vollständige Dokumentation der Programmausführung. 
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PROGRAMME ALS UNTERPROGRAMME 


Die Programme des Mathematik-Pakets können als Unterprogramme von Be- 
nutzerprogrammen im Speicher des HP-41C aufgerufen werden, Die speziellen 
Einspringmarken sind im Anhang B beschrieben. 


UMSPEICHERN VON PROGRAMMEN AUS MODULEN 


Wünschen Sie eine schrittweise Kontrolle bei der Programmausführung 
(TRACE), Modifikationen oder die Aufzeichnung auf Magnetkarten von Pro- 
grammen dieses Anwender-Moduls, so muß das Programm zunächst in den 
Speicher des HP41C kopiert werden. Erläuterungen zu der HP41C COPY- 
Funktion finden sich im Bedienungshandbuch. Für die Ausführung eines Pro- 
gramms ist es nicht notwendig, dieses Programm zu kopieren. 


UNTERBRECHUNG DES PROGRAMMLAUFES 


Alle Programme wurden so ausgelegt, daß sie richtig arbeiten, wenn sie von 
Anfang bis Ende durchlaufen, ohne den Rechner abzuschalten. (Beachten Sie 
dabei bitte, daß sich der Rechner auch selbst abschalten kann.) Wird der HPAIC 
abgeschaltet, kann es notwendig sein, daß für eine richtige Wiederaufnahme 
der Programmausführung das Flag 21 (SF21) gesetzt werden muß. 


MARKEN 


Es wird darauf aufmerksam gemacht, daß die Verwendung von ALPHA-Marken 
in einem selbstgeschriebenen Programm dann zu Problemen führen kann, wenn 
(diese Marken identisch mit den Marken im Programm aus dem Anwender-Modul 
sind. 
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MATRIXALGEBRA 


Dieses Programm berechnet Determinante und Inverse für Matrizen mit einer 
Dimension von maximal 14 x 14, wobei die Lösungen für Gleichungssysteme 
mit bis zu 14 Unbekannten bestimmt werden. 


Es wird das GAUSSsche Eliminationsverfahren mit modifiziertem Pivot-Ver- 
fahren verwendet. Aus Platzgründen kann das Verfahren nicht vollständig dar- 
gestellt werden, aber in der einschlägigen Literatur (vgl. Literaturangaben) kann 
die fehlende Information eingeholt werden. 


Das GAUSSsche Eliminationsverfahren besteht aus einer Folge von Zeilen- 
operationen zweier Typen: Die expandierende Elimination und die Lösung. 
Während der Eliminationsphase wird die N xN Matrix A in eine obere Dreiecks- 
matrix U umgewandelt, wobei A als nichtsinguläre Matrix angenommen wird. 
Die Multiplikatoren, aus denen diese Matrix erstellt wird, bilden eine untere 
Dreiecksmatrix L, welche 1.0 in der Diagonale enthält. Wenn auf das Pivot- 
Verfahren, das durch eine Folge von Zeilenvertauschungen zu einer Verbesserung 
der Genauigkeit bei vielen Gleichungssystemen führt, verzichtet wird, gilt für die 
Matrizen die Beziehung U = LA. Am Ende diese Programmteils wird Flag 4 
gelöscht, damit ist auch die Matrix A in ihrer ursprünglichen Form zerstört. Die 
ursprünglichen Elemente von A; sind durch die Elemente von U (iSj) und 
L (i>j) ersetzt. Die Rückschätzung in diesem Programm verwendet die trans- 
formierten Matrizen U und L für die Determinante und Inverse von A und zur 
Lösung des Systems homogener Gleichungen. 


Gleichungen: (am Beispiel einer 5 x 5 Matrix) 


EsseiA= I|A, An Ay Ay As 


Die Determinante von A, Det A, ergibt sich nach der Transformation nach U 
als Produkt der Diagonalelemente: 


Det A = (-1)* U), Uy U U U5;, 
wobei k die Anzahl der Zeilen ist, die vertauscht wurden. 
Es sei C die Inverse von A, d. h. diejenige 5x5 Matrix, die AC=CAe=1 erfüllt. 
Für Igilt, daß I; = 1 für i=j und 0 andernfalls. 
C wird spaltenweise in der folgenden Form bestimmt: 


Es sei c.;der j-te Spaltenvektor von C, z. B.: 
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ce; = C; I,j=1,2,3,4,5 


c.; wird aus der Lösung der Gleichung 


Ac.; =L.,, mit [.; = (1 für i>j, o andernfalls) 


Beispielsweise wird c., durch die Lösung von 


Ac.| = 


oOOOoOoo— 


bestimmt. 


Ein System von 5 Gleichungen für 5 Unbekannte kann geschrieben werden als: 
AyıXkı + A12%2 + A133 + A1aXg + AısXs = Bı 
AzıXı + A22X2 + AyaX3 + AyXy + AasXs = Ba 
AzıXı + A32%2 + A33%3 + AzyXy + AysXs = B; 
AyıXı + Ay2%2 + Ay3X3 + AyaXy + AysXs = By 
AsıXı + As2X2 + As3X3 + AsaXy + AssXs = B; 


wobei die x; unbekannt und die B; Konstante sind. 


In Matrixnotation wird diese Gleichung zu Ax = B, wobei X und B Spalten- 
vektoren sind: 


x] B, 
X) B, 
x= |%| undB= | B, 
X B, 
Xs B; 


Das Problem wird mit Ux = LB gelöst (wobei das Pivotieren vernachlässigt ist). 


Anmerkung: 


Ein Stop während der Programmausführung mit der Anzeige NO SOLUTION 
bedeutet, daß die Matrix A singulär ist. 
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@ Das Programm ist für 14x14 Matrizen ausgelegt, allerdings werden ab 7x7 
Matrizen Speichererweiterungs-Module benötigt. Konsultieren Sie die Matrix- 
Register-Karte für den Speicherbedarf verschiedener Matrizen. 


© Das Mathematik-Paket muß in eine Buchse nach dem Speichererweiterungs- 
Modul eingesetzt sein. 


® Sind die Elemente von A bereits richtig abgelegt, die Dimensionierung im 
Register 14 (R,,) enthalten, Flag 04 gesetzt und die Flags 06-10 gelöscht, 
dann kann die initialisierende Antwort ausgelassen und das Pivot-Verfahren 
durch Drücken von PVT gestartet werden. Soll ein Gleichungssystem 
gelöst werden, müssen Spaltenregister belegt und Flag 05 gesetzt sein. 


@® Wird DET als Unterprogramm aufgerufen, wird der Wert der Determinante 
ins Y-Register übergeben. 


@ BeiderEingabe der Matrixelemente darfdas Y-Registernichtgeändert werden. 


@® Die besten Ergebnisse werden erzielt, wenn die Matrix die Verfahrensvoraus- 
setzungen erfüllt. 


Literatur: 


Carnahan, Luther and Wilkes, Applied Numerical Methods, John Wiley and Sons, 
1969. 

George E. Forsythe, Michael A. Malcolm,and Cleve B. Moler, ComputerMethods 
in Mathematical Computation, Computer Science Department, Stanford Univer- 
sity, 1972. 

G.Forsythe and C. Moler, ComputerSolution of Linear Algebraic Systems, Prentice- 
Hall, 1967. 

C. Moler, “Matrix Computations with Fortran and Paging”, Comm. ACM, vol. 15, 
no. 4, pp. 268-270 (April, 1972). 


mr man Tre Tran Tanie 


Programm initialisieren [xEQ] MATRIX [ORDER=? 


Spalten/Zeilenzahl (N) eingeben R/S SET SIZE nnn 
(N=14) 


Umfang eingeben und weiter [xEQ] SIZEnnn 
[R/S] 


Matrixelemente zeilenweise eingeben 
(Ayı, Arz, Aız, .. Agı, Ayz Any .. [R/S] 
usw. [RS] 
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Schritt 4 wiederholen bis alle Elemente 
eingegeben 


Nach Bedarf: Matrix ausgeben 


Nach Bedarf: Matrix abändern ROW tCOL=? 


Eingabe von Zeile (I) und Spalte (J) 
des abzuändernden Elementes 
AlJ=? 


neuen Wert eingeben j AlJ= 
ROW# COL=? 


Schritte Nr. 7a und 7b wiederholen, 
wie benötigt. 


Ende Änderung 
Determinante berechnen 
Inverse bestimmen (spaltenweise) 


Für die Lösung simultaner Gleichungen, 
Spaltenvektor B eingeben 


Gleichungssystem lösen 


Nach Bedarf: Ausgabe der Spalte 


*Ein zusätzliches Drücken von wird benötigt, wenn kein Drucker angeschlossen ist. 
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Beispiel 1: 
Gesucht sind Determinante und Inverse der folgenden Matrix: 


ee | 
19.93 42 
23 4 2.8 
4 3.0.2, 
in 
Tasten: | Anzeige: 
MATRIX ORDER=? | 
5 (R7S) SET SIZE 50 (bedeutet, daß weniger als 


50 Register 
aktuell verfügbar sind) 


SIZE (ALPHA) 050 

A1,1=? 
6[R/S) 3[R/S) 2[CHS) [R/S) 2[R7S) 

3[R7S) 1[R/S) 4[R/S) 3[cHs) [R7S) 

4[R/S) 2[R/S) 2{R7S) 3 [R/S) 

1 (cHS) [R/S) 2[cHs) [R/S) 9 (RZ) 

4{R/S) 3{R7S) O[R/S) 2[R7S) 1[R7S) 

3[R/S) 5[{R7S) 6[CHs) R/S) 


6 2 0.0000 

DET DET=-200.0000 (Det A) 
INV C1, 1=0.2000 
C2, 1=-1.7500 
C3, 1=0.7000 
C4, 1=1.7250 
C5, 1=1.0000 
R/S C1, 2=-0.1200 
C2, 2=-2.0000 
R/S C3, 2=1.2800 
C4, 2=2.5400 
R/S C5, 2=1.4000 
C1, 3=-0.0400 
C2, 3=0.5000 
C3, 3=-0.2400 
C4, 3=-0.5700 
R/S C5, 3=-0.2000 


2 
Den 
7) 


e.] 
_ 
17) 


a 12] |2] 121 12921120 
<= I=1 =] I>=] In 
ol Io} In Io] Io] Io 


bi) 
— 
9 


Matrixalgebra 
Anzeige: 


C1, 4=-6.0000 E-10 
C2, 4=1.7500 

C3, 4=-0.5000 

C4, 4=-1.6250 

C5, 4=-1.0000 

C1, 5=-1.0000 E-10 
C2, 5=1.5000 

C3, 5=-1.0000 

C4, 5=-1.7500 

C5, 5=-1.0000 


15 


Unter Anwendung der Schleifentechnik auf den untenstehenden Schaltkreis sind 
die Ströme ],,I,, I; und I, zu bestimmen. 


10 109 19 
NN — MN — AN 
ei A pe 
34V - 210 10 zn sin 
| 1 3 4 
ine BEER \ 
Die zu lösenden Gleichungen lauten: 
21, -I = 34 
-], +3], Tu 0 
-J +3]; -J = 0 
+3], = 0 
In Matrixnotation: 
2-00 I; 34 
1 3 1 0 9 0 
0- 3-1 I; 0 
00-3 Lu 0 


16 _Matrixalgebra 


Tasten: 


MATRIX 
4 (R/S) 


2[R/S) 1[c#s) [R7S) O[R/S) O[RS) 
1 [EHs) [R7S) 3 [R7S) 1 (cHs) [R7S) 
o[R/S) o{R7S) 1[cHs) [R/S) 3[R7S) 
1 [CHS) [R/S) 0 [R/S) 0 [R7S) 

1 (EHs) (R7S) 3 [R/S) 

KXEQ) (ALPHA) SIMEQ [Arena] 

34 (R/S) O[R/S) O[R/S) 0 (R7S) 
R/S 


Anzeige: 


ORDER=? 


A1,1=? (es wird weiterhin ein 


Umfang von 050 
angenommen.) 


0.0000 
B1=? 
0.0000 
X1=21.0000 
X2=8.0000 
X3=3.0000 
X4= 1.0000 


17 
LÖSUNGEN FÜR F(x)=0 IN EINEM INTERVALL 


Dieses Programm verwendet einen modifizierten Sekanten-Algorithmus, um 
reelle Lösungen der Gleichung f(x)=0 zu bestimmen. Der Anwender muß die 
Funktion eingeben und kann zwei Startwerte (x, und x,)angeben, die nahe an der 
gewünschten Lösung liegen sollen. Wird das Intervall anfangs nicht definiert, so 
wird ein Anfangsintervall zwischen | und 10 angenommen. 


Gilt f(x))-f(x,) = 0 und ist die Funktion im Intervall stetig, so wird immer eine 
Lösung gefunden. Gilt dagegen f(x,)-f(x>) > 0, kann die Suche nach einer Lösung 
erfolglos sein. Existiert mehr als eine Lösung im betrachteten Intervall, so wird 
eine Lösung bestimmt und der Benutzer kann ein kleineres Intervall definieren 
und das Programm wiederholen. 


Die Funktion f(x) kann unter jeder beliebigen globalen Marke - mit höchstens 
6 Zeichen Länge - eingegeben werden, es wird jeweils x zu Anfang im X-Register 
erwartet. Es können mehrere Funktionen in den Programmspeicher geladen 
werden,dader Name der zu bestimmenden Funktion vom Programm erfragt wird. 
Das Programm belegt die Register 00-06. Die restlichen Register und der Stack 
sind für die Definition von f(x) frei. 


18 Lösungen für f(x) in einem Intervall 


Umfang: 007 


| ümmer | Anweisung —__Teingete [runtn _Ianeie 


Funktionseingabe vorbereiten 


[PRGm ]-Mode einstellen, Funktion 
unter dem gewünschten Label ablegen; 


hinzufügen und [[PRGM_]-Mode 
ausschalten. 


Prögramm initialisieren FUNCTION NAME? 


Funktionsnamen eingeben, (ALPHA- 
Modus ist in Schritt 3 gesetzt) 


Falls zwei Anfangsschätzungen einge- 
geben werden sollen, ersten Wert ein- 
geben, ansonsten weiter mit Schritt 7 


zweite Anfangsschätzung eingeben 
Programm ausführen 


um eine andere Lösung zu bestimmen ist 
zu drücken. Für weitere Be- 
rechnungen weiter mit Schritt 1 oder 3. 
Anm.: Es gibt drei Arten der normalen 
Beendigung des Programmes, Es 
erscheinen dabei die folgenden 
Meldungen: 
1) keine Lösung gefunden 
2) Lösung ist <Wert> 
3) Lösung liegt zwischen <Wert 1> 

und <Wert 2> 
Bei Verwendung als Subroutine (sofern 
die Anfangsschätzungen bereits ein- 
gegeben), können die Schritte 3-9 durch 
Drücken von SOL übersprungen 
werden. 


GUESS 1=? 


GUESS 2=? 


(vgl. Anm.) 
(vol. Anm.) 


Lösungen für f(x) in einem Intervall 19 
Beispiel 1: 


Gesucht sind die Lösungen für In x + 3x - 10.8074 = 0. f(x) ist durch FF 
definiert. 


Tasten: Anzeige: 


SIZE 007 
Gojaleı 


Fr AR) 
KJEIIEZ 

10.8074) 

RTN 


PRCSA 

(ALPHA) SOLVE FUNCTION NAME? 
FF GUESS 1 =? 

R7S ROOT IS 3.2134 


Es ist derjenige Winkel « zwischen 100 und 101 Rad gesucht, für den sin « = 
0.01 ist. Somit ist f(x) = sinx - 0.01. Verwenden Sie dabei [LBLJ ANGLE. 


Tasten: Anzeige: 


Galalcı 
ANGLE 
RAD 


[sın) 

Klcı 

DEG 

RTN 

SOLVE FUNCTION NAME? 

ANGLE GUESS 1 =? 

100 GUESS 2 =? 

101 ROOT IS BETWEEN 100.5410 
R/S AND 100.5410 


u mehr signifikante Ziffern angezeigt werden, drücken Sie [FX] 9 und 


20 Lösungen für f(x) in einem Intervall 
Beispiel 3: 
Gesucht sind die Lösungen von x? + 1 = O unter Verwendung von CC 


Tasten: Anzeige: 
GC) 
66 
&2 ı 

RTN 
SOLVE FUNCTION NAME? 
cc GUESS 1 =? 

R/S NO ROOT FOUND 


21 
NULLSTELLEN UND FUNKTIONSWERTE VON POLYNOMEN 


Dieses Programm kann zur Bestimmung von reellen Nullstellen bei Polynomen 
höchstens fünften Grades verwendet werden, wenn der Koeffizient des Terms 
mit höchstem Grade 1 ist. Die Gleichung lautet demnach: 


rt Farrell n=2,3,4or5. 


Ist der führende Koeffizient ungleich 1, so muß - durch Division der gesamten 
Gleichung mit diesem Koeffizienten - dieser zu 1 gemacht werden. 


Beim Initialisieren des Programms ist der Polynomgrad (n) vom Benutzer anzu- 
geben. Der Rechner erfragt dann die Koeffizienten a,., ... , A, &. Sind die 
Koeffizienten gleich Null, so ist dies anzugeben. Die Koeffizienten sind in den 
Registern 00-04 gespeichert. 


Gleichungen: 


Die Routinen für Gleichungen vierten und fünften Grades verwenden ein 
Iterationsverfahren, um eine reelle Lösung der Gleichung zu finden. Für diesen 
Fall darf a, nur ungleich 0 sein, (Ist a,= 0, dann ist Null eine reelle Lösung und da- 
mit kann der Grad des Polynoms durch Abdivision von x um 1 reduziert werden.) 
Sobald eine Nullstelle gefunden ist, wird die Gleichung durch synthetische Divi- 
sion auf eine Gleichung zweiten oder vierten Grades reduziert. 


Zur Lösung der Gleichung vierten Grades muß zuerst die kubische Gleichung 


»rbyrbytrb=0 


gelöst werden, wobei 1, = -% 
b, = a3aı = 4%, 
bo, = % (4a, — a2) — ar ist. 


Es sei y, die größte Lösung der obigen Gleichung dritten Grades. 


Damit kann die Gleichung vierten Grades in zwei quadratische Gleichungen zer- 
legt werden: 


*+(A+Cx+(B+D)=0 

"+ A-x+B-D,=0 
wobei A=—, -2 ‚D=VB-„.,C= VA-3,+y 
ist. 


22 Nullstellen und Funktionswerte von Polynomen 


Lösungen für die Gleichung vierten Grades ergeben sich aus der Lösung der 
beiden quadratischen Gleichungen. 


Eine quadratische Gleichung x? + a‚x + a, = 0 wird nach der Formel 


Mc ar öst. Ist D = ar _ > 0,sind die Lösun- 
_ >: 2 ai a. gelöst. Ist 7 a) „Ss 


gen reell; ist D<O, sind die Lösungen komplex nach: u + iv= — — =ZiV-D. 


Eine reelle Lösung wird als einzelne Zahl ausgegeben. Komplexe Nullstellen er- 
scheinen immer als Paare der Form utiv und sind in der Ausgabe bezeichnet. 


Anmerkung: 


® Lange Programmlaufzeiten sind für Gleichungen 3.-, 4.- oder 5-ten Grades zu 
erwarten, da diese ein- oder mehrfach eine Iteration durchlaufen. 


® Das Programm belegt die Register 00-22 


Programm initialisieren DEGREE=? 


Grad des Polynoms eingeben a(n-1)=? 
(n=2, 3, 4, 5) 

Koeffizienten a,-1 eingeben 
Koeffizienten= 0 müssen als 0 ein- 
gegeben werden, 

Schritt wiederholen, bis a9 

In der Anzeige erfragt wird, 


% (n-2)=? 


Koeffizienten a, eingeben 


Für Lösungen ist [R/S. zu drücken, es 
werden dann alle Lösungen richtig be- 
zeichnet. Weiter mit Schritt 9 


für den Funktionswert ist mit Nein (N) zu 
antworten 


x eingeben, f(x) wird angezeigt 
für ein neues x; x eingeben 
und[RYS Jdrücken 


für ein neues Polynom gleichen Grades 
weiter mit Schritt 1 oder andern der 
entsprechenden Koeffizienten in den 
Registern 00-04 und ROOTS 


Nulistellen und Funktionswerte von Polynomen 


MBaionscl . 
beispiei 1: 


Gesucht sind die Nullstellen von x° - x*- 101x? + 101x? + 100x- 100 =0 


Tasten: 


SIZE (ALPHA) 023 
POLY 

5 [R/S) 

1 [EHs) [R7S) 

101 (CHs) [R7S) 

101 [R/S] 

100 (R7S) 

100 [C#s) [R7S) 


R/S 


Anzeige: 

DEGREE=? 

a4=? 

a3=? 

a2=? 

al=? 

a0=? 

ROOTS? 

ROOT=10.0000 (Lösung 1) 
ROOT= 1.0000 (Lösung 2) 
ROOT=1.0000 (Lösung 3) 
ROOT=-1.0000 (Lösung 4) 
ROOT=-10.0000 (Lösung 5) 


Es ist 4x? - 8x? - 13x? - 10x + 22 = 0 zu lösen. Die Gleichung muß zu 


4 » BB. 1 22 
- 2x - —y’- —y + 
x X A X A X 4 
Tasten: 
[XEa)[ALPHA)poLYy (ALPHA) 
4 RB) 
2 (C#3) RS) 


13 Eure) 4 G)[cHs) (RS) 
10 Eee) 4[+)[cHs) (R7S) 
22 [ENTeRe) 4 (+) [R7S) 


= () umgeformt werden. 


Anzeige: 


DEGREE=? 
a3=? 
a2=? 
al=? 
a0=? 
ROOTS? 
=-1.0000 
V=1.0000 
=-1.0000 


sind -1.00 + 1.00 i) 


-V=-1.0000 


ROOT=3.1180 
ROOT=0.8820 


(Lösung 3) 
(Lösung 4) 


23 


(Die Lösungen 1 und 2 


24 Nullstellen und Funktionswerte von Polynomen 


Beispiel 3: 
Wie verändern sich die Lösungen, wenn der Koeffizient von x? von -13/4 in -5 
geändert wird? 


Tasten: Anzeige: 
5[cHs) 02 . 
ROOTS U=-1.1386 (Lösungen 1 und 2 sind 
U=-1.1386 
-V=-0.8555 
R/S ROOT=3.5031 (Lösung 3) 
R/S ROOT=0.7741 (Lösung 4) 
Beispiel 4: 


Es sind die Funktionswerte des Polynoms 


f(x) = x° + 5x’ - 3x? - 7x + 11 an den Stellen x=2.5 und x=-5.0 zu bestimmen. 


Tasten: Anzeige: 

POLY DEGREE=? 

5 a4=? 

5 a3=? 

0 a2=? 

3[CHs) al=? 

z[cHs) a0=? 

11 ROOTS? 

N X=? 

2.5 F<X>=267.7188 


5[EHs) F<X>=-29.0000 


25 


NUMERISCHE INTEGRATION 


- 
x; ) -TRAP +SIMP 


Dieses Programm leistet numerische Integration für den Fall, daß eine Funktion 
explizit oder nur an einer endlichen Zahl abstandsgleicher Stützstellen (diskreter 
Fall) definiert ist. Die Integrale expliziter Funktion werden nach der 
SIMPSONschen Regel bestimmt; für diskrete Fälle wird entweder nach der 
Trapezregel angenähert oer auch die SIMPSONsche Regel angewendet. 


Diskreter Fall 


Es seien X, X X34 =. , Xn n abstandsgleiche Punkte (mit xj= xo+jhj=12 ... ,n), 
denen dir bekannten Funktionswerte f(x,),[(X 1), -.- , [(%,) zugehören. Die Funk- 
tion f(x) selbst muß nicht bekannt sein. Nach der Eingabe der Schrittweite h 
und der Funktionswerte ((x,)j=0, 1,2, ...,n wird das Integral 


f dx 
% 


entweder nach der Trapezregel angenähert: 
Xn h n-1 
f f(x) dx | +2 > f(x) + ) 
Xo jel 


oder nach SIMPSON geschätzt: 
Kn 
Ix)dı = aß + 4f(x,) + 2f(x,) + 
% 
+ 4fl&n-3) + 2fin-2) + 4n-ı) + | 


Um die SIMPSONsche Regelanwenden zu können,muß ngerade sein. Istn unge- 
rade,stoppt der Rechnermitder AnzeigeN NOT EVEN, sobald [D] gedrückt wird. 


26 Numerische Integration 


Explizite Funktionen 


Ist die Gestalt der Funktion f(x) bekannt, so kann die Funktion im Programm- 
speicher abgelegt und nach SIMPSON numerisch integriert werden. Der Benutzer 
muß die Intervallgrenzen a und b, sowie die Zahl der Teilintervalle n ‚in die (a,b) 
zerlegt wird, angeben. Dieses n muß gerade sein, andernfalls wird NNOT EVEN 
angezeigt. Das Programm wird fortgesetzt mit 


o=-3,%,=-x%+jh,j=1,2,..,n- l,und x, = b, wobei 


Das Integral f f(x)dx wird nach der obigen SIMPSONschen Regel 
angenähert. a 


Die Funktion f(x) kann unter Verwendung einer beliebigen, maximal 6 Zeichen 
langen, globalen Marke in den Programmspeicher eingegeben werden, wobeix im 
X-Register erwartet wird. Es können mehrere Funktionen in den Programm- 
speicher eingegeben werden, da der Funktionsname vom Programm vorher er- 
fragt wird. Das Programm belegt die Register 00-07; die übrigen Register stehen 
für die Definition von f;(x) zur Verfügung. 


Umfang: 008 


Maske auf Rechner legen 


für explizite Funktionen weiter 
mit Schritt 9, für den diskreten 
Fall weiter mit Schritt 3 


DISKRETER FALL 
Programm initialisieren 


Abstand zwischen den x-Werten 
eingeben 


Funktionswerte für x; eingeben 
Schritt für j=0,1, ... n wieder- 
holen 


Fläche nach Trapezregel bestimmen 


Fläche nach SIMPSON bestimmen 
(n gerade) 


für weitere Berechnungen weiter 
mit Schritt 2 


EXPLIZITE FUNKTIONEN 
Funktionseingabe vorbereiten 


KXEOJINTG 


A] 


D 


Numerische Integration 27 


nn 


In [LPRGM_]-Modus schalten; Funktion 
unter gewünschter Marke eingeben; nach 
[RTN)wieder PRGM-Modus ausschalten. 


Programm initialisieren 
Intervallgrenzen a und b eingeben 


Anzahl der Teilintervalle n eingeben 
(n gerade) und Fläche nach SIMPSON 
berechnen. 


Funktionsnamen eingeben 

sollen a, b oder n geändert werden, mit 
entsprechenden Schritt fortsetzen; für 
weitere Berechnungen nach Schritt 2 
gehen. 


ı ik. 


Gegeben seien die untenstehenden acht Werte für f(x), j= 0,12,..,8; gesucht 
sind die Näherungen des Integrals 


f j x) dx 


nach Trapez- und SIMPSONscher Regel. Der Wert für h ist 0.25. 


35 Hoss, ELIEm 
EmE ; 2.8 73 . 7 57 71116 


Tasten: Anzeige: 


SIZE 008 
INTG 0.0000 


28 Numerische Integration 


Tasten: 


.25[4) 2(B] 
2.8[B) 3.8(®) 
5.2(6) 7(e) 
9.2(B) 12.1(E) 
15.6(B) 20(8) 
[e) 

(D) 


Beispiel 2: 


Gesucht ist der Funktionswert von 


Anzeige: 

16.6750 (Trapezoidal) 

16.5833 (Simpson’s) 
dx 


1-cosx + 0.25 


für n=10 und n=30. Man beachte, daß x in Rad erwartet wird. Es ist eine gute 
Gewohnheit, am Anfang der Routine den Winkelmodus auf Rad und am Ende 
der Routine wieder auf Grad zu schalten, wenn hauptsächlich in Grad gearbei- 
tet wird. Die Funktion soll unter FF eingegeben werden. 


Tasten: 


GOlKlER) 

FF (ALPHA) 

RAD [ALPHA ) 
f 

EIBZIEI 

256) 17%) 

DEG 
RTN 

(xEa) (ALPHA) INTG [ALPHA 
VEL=SOWIEZEIGIEN 

10mE® 

FF (R7S) 

som (le) 

FF (R7S] 


Anzeige: 


FUNCTION NAME? 
8.2193 (n=10) 
FUNCTION NAME? 

8.3774 (n= 30) 


Die exakte Lösung ist Eu = 8.3776 


29 
DIFFERENTIALGLEICHUNGEN 


Dieses Programm berechnet Ableitungen ersten und zweiten Grades nach der 
Methode von RUNGE-KUTTA. Eine Ableitung ersten Grades ist von der Form 
y’=f(x,y), mit den Anfangswerten x,,y, und y,. 

In jedem Fall kann die Funktion f(x) in den Programmspeicher unter einer 
beliebigen, maximal 6 Zeichen langen, globalen Marke eingegeben werden, 
Die Werte x und y werden im X- bzw. Y-Register erwartet; y’ wird im Z-Register 
bei Ableitungen zweiten Grades erwartet. Das Modul-Programm belegt die Re- 
gister 00-07. Die verbleibenden Register stehen für die Definition der Funktion 
zur Verfügung. 


Die Lösung ist ein numerischer Wert, beidemy; fürx;=x,+tih(i=1,2,...) fürein 
vom Benutzer zu bestimmendes Inkrement berechnet wird. Der Wert für h kann 
zu jeder Zeit während der Programmausführung geändert werden, man speichert 
dazu h/2 in das Register 01. Dies erlaubt Lösungen der Gleichung beliebig nahe 
an einem ihrer Pole (y> +»). 


Gleichungen: 


Erste Ableitung: 


ud + + 20, + 2c3 + C,) 


wobei 
c, = hf(x,, yı) 
c 
= (x + —, yı + — 
C2 
G = (x +, yı 5 
s=hfx, +h,yı + ca) 
Zweite Ableitung: 


yaznsn| + + k, + «) 


Yin yi + r (kı + 2k2 + 2ks + ku) 


30  Differentialgleichungen 
kı = hf &, Yu, yı) 


h E,.& PR: 
ka us 7 7’ yı Br yı + gKn RE T 


h Rn Fa 
ka (x T zZ yı Wr Ya gkı yı 2) 


| A n(« + h, yı + hy; + > ka, yı + «) 


Anmerkung: 


@ Werden die Werte für eine zweite Ableitung eingegeben, so müssen die 
Werte von x, und y, vor dem Wert für y, eingegeben werden. Alle Werte 
müssen eingegeben werden, auch wenn sie gleich Null sind. 


Umfang: 008 


Funktionseingabe vorbereiten f(x,y,y') 


PRGM-Modus einschalten, Funktion 
unter der gewünschten Marke ein- 
geben; nach den PRGM-Modus 
ausschalten. 


Programm intialisieren DIFEQ FUNCTION NAME? 
Namen der Funktion eingeben Name ORDER=? 
Eingabe ob 1. oder 2. Ableitung 1 oder 2 STEP SIZE=? 

(1 oder 2) 

Schrittweite h eingeben X0=? 
Anfangswert für x eingeben Yo=? 
Anfangswert für y eingeben 


für eine zweite Ableitung 
(Abl. 2. Ordnung) Anfangswert für y' 
eingeben. 


Ausgabe der Werte für x und y 
nacheinander 


Differentialgleichungen 31 
Beispiel 1: 
Unter Verwendung von EX, soll die Differentialgleichung erster Ordnung 
‚_sinx+ tan”! (y/x) 
a In(Vx? + y?) 
numerisch bestimmt werden, wobei x,=y,=1. Der Winkelmodus muß auf RAD 
gesetzt werden. Außerdem werden drei weitere Speicherregister für die Definition 
der Funktion benötigt. 
Tasten: Anzeige: 


SIZE (ALPHA) 011 
ro)LIL) 

[PRGM) 

Fx [atHa) 

eo) (ALPHA) RAD (ALPHA) 

08 


oO 
© 


zlefzielejelglälel: 
I vi I] Jar 
oO —_ 
[0 +) oO 


2 
oO 
r 
—- oO 
te) 


SUR 
oO 


DEG 


RTN 


[xEa) [ALPHA )DIFEQ FUNCTION NAME? 

FX [R7S) ORDER=? 

1 STEP SIZE=? 

‚5{R/S) X0=? 

1 Yo=? 

1 1.5000 x.) 
2.0553 (yı) 
2.0000 (x) 


32  Differentialgleichungen 


Tasten: Anzeige: 

R/S 2.7780 (y2) 
2.5000 (X3) 
3.2781 (y3) 

etc. 

Beispiel 2: 


Unter Verwendung von DIF, soll die Gleichung 2. Ordnung 
d- Wi Ta 


gelöst werden, wobei X9=yo=yo=0 und h = 0.1. 


Man formuliert die Gleichung um zu: 


un U = Dre) xı#l 


1- x | 


Tasten: Anzeige: 


Goleslee 

DIF [ALPHA ) 
08 

1E) 

08 

E3) 

&? 
IE} 


DIFEQ FUNCTION NAME? 

DIF ORDER=? 

2[R’s) STEP SIZE=? 

.1 X0=? 

0 Yo=? 

0 Yo.=? 

0 0.1000 (x)) 


Anzeige: 


0.0002 
0.2000 
0.0013 
0.3000 
0.0046 
0.4000 
0.0109 
etc. 


Differentialgleichungen 


(yı) 
(X2) 
(y2) 
(X3) 
(Ys) 
(%4) 
(y4) 


34 
FOURIER-ANALYSE 


Jede periodische Funktion kann - unter Verwendung der nachfolgenden Formeln 
- als Folge von Sinus- und Kosinusfunktionen dargestellt werden: 


fit) = >- + 3 (ve cos Br b, sin 2x) 


kl 


— > + TG cos E= ) 


k=1 


7 
a, = 2[ f(t) cos Zutk d,k=0,1,2, 
T% T 


= an ft) sin FE d1,x= 1,2,. 


= (ag + by)” 


0,= RS 
Ak 


T = Periodizität von fft) 


Dieses Programm bestimmt die Fourier-Koeffizienten für diskrete Fälle der 
obigen Formeln, wenn eine große Zahl von Stichproben der periodischen Funk- 
tion gegeben ist. Bis zu zehn aufeinanderfolgende K.oeffizientenpaare könnenauf 
einmal aus N abstandsgleichen Punkten berechnet werden. Die Koeffizienten 
können sowohl in kartesischen Koordinaten als auch in Polarform dargestellt 
werden. 


Der Wert für N sollte so gewählt werden, daß er mindestens doppelt so groß ist wie 
der höchste erwartete Wert der Grundfunktion, die in der periodischen Funktion 
vorliegt. Eine zu geringe Schätzung für N bewirkt, daß eine Hälfte des Stich- 
probenintervalls mit einer niedrigeren Frequenz überschätzt wird (ein Umstand 
der als ’aliasing’ bekannt ist). 


Die Register 00-26 sind vom Programm belegt. 


Programm initialisieren 

Anzahl der Stichproben in einer 

Periode eingeben #Stichp. 
Anzahl der gewünschten Frequenzen 

eingeben # Frequ. 
Exponent der ersten Funktion (J) 

eingeben J 

Yn mit n=1, 2, ... , N eingeben Yn 


Schritt 5 wiederholen bis RECT? 
in der Anzeige 


Ist die Antwort ja, muß gedrückt 
werden, um die Koeffizienten für 
JSksJ+ #Frequ. in kartesischen 
Koordinaten auszugeben. 


Ist die Antwort Nein (N), sollen N 

die Koeffizienten in Polarkoordinaten 
ausgegeben werden. 

durch drücken von werden auf- 
einanderfolgende Koeffizienten angezeigt, 
um den Wert der Fourierreihe zum 
Zeitpunkt t zu bestimmen, ist der 
USER-Modus zu setzen und t 
einzugeben. 


Fourier-Analyse 


NO.FREO=? 
1ST COEFF=? 


Yie..RECT? 
Y2=?..RECT? 
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Beispiel: 

Es ist eine Fourier-Serie für die unten gezeigten Wellenformen zu bestimmen. Da 
12 Stichproben vorliegen, werden 7 Frequenzen ausgewählt (DC-Term und 6 
harmonische Funktionen). Die Koeffizienten sind als Polarkoordinaten darzu- 


stellen. 


20 
15 
10 
5 
0 
-5 
-10 
-15 
t £(t) 
l 14.758 
2 17.732 
3 2 
4 «12, 
5 - 7.758 
6 -11 
7 - 9,026 
8 12, 
9 2 
10 14.268 
11 10.026 
12 15 
Tasten: Anzeige: 
SIZE 027 
FOUR (ALPHA) NO. SAMPLES-? 
12 [R/S) NO. FREQ=? 
7{R/S) 1ST COEFF=? 


O[R/S) Y1=? 


Tasten: 


14.758 [R/S) 
17.732 [R/S) 

2 [R/S) 

12 (CHs) [R7S) 
7.758[CHs) (R/S) 
11[cHs) [R/S) 
9.026[CHs) (R’S) 
12 [Chs) [R/S) 
2[R’S) 

14.268 [R/S) 
10.026 [R/S) 

15 (R/S) 


R/S 


2 
— 
2) 


2 
— 
172) 


2 
— 
1) 


2 
= 
[77] 


2/21 |2) |2/} /|2/|2 
SI ISIS IS IS In 
o||Iol io] Io] In] io 


2| /|2/|2 
I IS IS 
o| Io Io 


Somit ist flt)= 2+ 15 cos 


Art 


Er. 
sın 12 


t 
5 cos 72 


1Ort 


+ 3:c08 


Zrrt 


+ sin 


Fourier-Analyse 


Anzeige: 


Y2=? 

Y3=? 

Y4=? 

Y5=? 

Y6=? 

Y’=? 

Y8=? 

Y9=? 

Y10=? 

Y11=? 

Y12=? 

RECT? 
a0=4.0000 
b0=0.0000 
a1=14.9998 
b1=1.0000 
a2=3.0000E-8 
b2=1.0000 
a3=-5.0000 
b3=1.0000 
a4=3.3333E-9 
b4=3.2000E-9 
a5=3.0002 
b5=1.4673E-5 
a6=0.0000 
b6=2.3599E-8 


Zrrt 
12 


Ort 


m 
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KOMPLEXE OPERATIONEN 


Dieses Programm gilt für verkettete Operationen mit komplexen Zahlen in kar- 
tesischer Form. Neben den vier Grundrechnungsarten fürkomplexe Zahlen (+, -, 
x,:) sind die gebräuchlichsten Funktionen der komplexen Zahlen z und w (|zl, 
l/z, z", z’" e” 1nz, sinz, cosz, tanz, a”, loga*, z”” und z”) enthalten. Funktionen 
und Operationen können gemischt werden, um Ausdrücke wie 23/(zı+23), 
e22 iz, + z,| +|z, - z,|etc., zu berechnen, wobei z|, z, und z; komplexe 
Zahlen der Form x+iy sind. 


Für den wiederholten Gebrauch dieser Operationen, kann der Benutzer die ein- 
zelnen Programme ausgewählten Tasten zuweisen und eine entsprechende 
Maske selbst herstellen. Eine vernünftige Neuzuweisung könnte etwa sein: 


(asn)sınz [sın) 
[asn) LNZ 
(asn) C+ 
Gsn)c- 2) 
(Asn) CıNV 


Das Logik-System für diese Programme kann man sich am besten als eine Art 
umgekehrte Polnische Notation mit einem aus zwei komplexen Zahlen be- 
stehenden Stack vorstellen. Das untere und das obere Register des komplexen 
Stacks sei als & bzw. r bezeichnet. Diese entsprechen somit den X- und T-Re- 
gistern des Vier-Register Stacks des Rechners*. Eine komplexe Zahl z; wird 
durch die Tasten y, (ENTER?) x, ins E-Register transportiert. Durch die Eingabe 
einer zweiten komplexen Zahl z, ([ENTERt] y, x,) wird z, ins Register ı 
geschoben und z, im Register E abgelegt. Der vorherige Inhalt von t ist verloren. 


Funktionen bedienen sich des E-Registers, deren Ergebnis im E-Register steht (mit 
Ausnahme von]|z|, wobei eine reelle Zahl erscheint). Arithmetische Operationen 
beziehen sich auf das &- und r-Register; das Ergebnis wird in & abgelegt. 


Das Programm des Anwender-Moduls belegt die Speicher 00-04. 
Gleichungen: 


Es sei 


4=%xtiy=nef,k=1,2 
z=x+iy=re® 
* Jedes Register des Stacks muß zwei Zahlen enthalten: den reellen und den imaginären Teil des 


komplexen Inhalts. Somit werden zwei Stack-Register für eine komplexe Zahl gebraucht. Komplexe 
Stack-Register werden (sprachlich) wie ein einziges Register behandelt. 


Komplexe Operationen 
Die Lösung sei für jeden Fall von der Form u + iv. 


z + 2=(X +x)+ ilyı + ya) 


Zı = 22 = (Xı = %e) + i(yı — yo) 


Z, 2,7 rırge1: + 92) 
Zul 2% = u ei, - 02) 
IR 
|z| = Vr®+ y2 
l/z N i I 
r re 
zR = rneind 
(2+ ze) 
zn = line \n 2 7/ k=0,1l,.,0r1 


(es werden alle Lösungen ausgegeben, k=0, 1, ... n-1) 
e?= e* (cos y + isin y), mity in Rad 
Inz = Inr + i6, mit z#0 


a? = ezina mit a> 0 undreell 


log.z = Ann mit a>0 reell und, z# 0 
Ina 


zW = ewinz mit z# 0, w is komplex 
zw = elnziw mit z# 0, und komplex w # 0 
sinz = sinxcoshy + icosxsinhy, Winkel in Rad 
cosz = cosxcoshy — isinxsinhy, Winkel in Rad 


sin 2x + sinh2y 


———, Wi in Rad 
cos 2x + cosh2y MADRE N 


tanz = 
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Umfang: 005 


ARITHMETIK 
erste komplexe Zahl eingeben 


(X +iy) 


ENTER+ 


zweite komplexe Zahl eingeben 
(X + iyp) 


auswählen einer der vier Operationen: 
@ Addieren (+) 


@ Subtrahieren (-) 


® Multiplizieren (x) 


@ Dividieren (*) 


das Ergebnis der Operation ist nun im 
Stack; weiter mit Schritt 5 für Funktionen 
oder Schritt 2 mit anderen arithmetischen 
Operationen 


FUNKTIONEN 
Auswahl einer der folgenden Funktionen: 
@ Absolutwert (z) 


x) 
m 
ie} 
— 
> 
N 


| 


@ Reziprokwert (1/z) 


mie 
m \ 
kejiesit 
_ ım 
<> ID 
=|+ 
= 


25 


® ganzzahlige Potenz von z 
4) 
KEOJZ N 
R/S 
© n-te Wurzel von z (z"’") 
Beachte: es gibt n Lösungen ENTER+ 
der Form u+iv ZrAIN 
eo. 
KXEQJer ZZ 
@ natürlicher Logarithmus von 
2 (Inz) LNZ 
@ z-te Potenz einer reellen 
Zahl a (a?) 
ENTER+ 


B 
m 
5 
{2°} 
nr 
N 


“ 
Ö 
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@ Logarithmus von z zur Basis 
a (l0g,”) 


@ Komplexe Potenz von z 
w=%o + iya (z“) 


© w-te komplexe Wurzel von z 


(er 


KXEQ]LO 
R/S 
KXEQ]7 
[RS] 
KXEQ]Z 


© sin(z) 


> 
N 


IE 


© cos(z) 


m 

(02) 
DD 
oO 
oo 
"N 


23 


@ tan(z) 


o 
2 


NZ 


weiter mit Schritt 2 für arithmetische 
oder Schritt 5 für andere Operationen 


Beispiel 1: 
Es soll der Ausdruck 
Be 
Ztz ' 
mit z,=23+13i, z,=-2+i und z3=4-3i berechnet werden. 


(Vorschlag: danurjeweils zwei Zahlen vom Programm bearbeitetwerdenkönnen, 
sollte die Gleichung umformuliert werden zu 


x [1/, + z3)].) 
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Tasten: 


SIZE [ALPHA] 005 


1 (ER) 2[cHs) [EnteRe) 
3 [cHs) [EmeRe) 4 

(XEo) [ALPHA] C+ (ALPHA) 

[xEa) [ALPHA] Cınv [ALPHA] 
R/S 

13(EER) 23 

cx [ALPHA] 


Beispiel 2: 


Anzeige: 


U=2.0000 
V=-2.0000 
U=0.2500 
V=0.2500 


U=2.500 
V=9.0000 


Gesucht sind die drei dritten Wurzeln von 8. 


Tasten: 


0 (EıER) 8 (Eımre) 3 
[XEQ) [ALPHA ) Zr 1/N (ALPHA) 


Beispiel 3: 


2 
Berechne e? ° ‚mit z=(1Hi): 


Anzeige: 


U=2.0000 
V=0.0000 
=-1.0000 
V=1.7321 
=-1.0000 
V=-1.7321 


Anzeige: 


U=0.0000 
V=2.0000 
U=0.0000 
V=-0.5000 
U=0.8776 
=-0.4794 


real (z, + z,) 
imag (2, + z,) 
1/ (Z, ns Z;) 


(zı/(z, + Z3)) 


(z?) 
(z°?) 
(e2?) 
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Beispiel 4: 

Berechne sin(2 + 3i): 

Tasten: Anzeige: 
3 (ER) 2 

SINZ [ALPHA U=9.1545 


=-4,1689 


HYPERBOLISCHE FUNKTIONEN 


Dieses Programm berechnet die hyperbolischen Funktionen und deren Inverse. 
Falls die einzelnen Programme bestimmten Tasten des Rechners zugeordnet 
werden sollen, um sie mit einer entsprechenden Maske zu kennzeichnen, könnte 
eine Umdefinition von Tasten etwa so aussehen: 


SINHESIN) 
COSHLCOS) 
TANHLTAN) 
ASINHBEIN) 
ACOSHILEOS) 
ATANHDCTAN) 


aalalaa]e 


Gleichungen: 
Hyperbolische Funktionen: 


mn [et] 
eN 


x -X 
coshx=-£ Eng 
2 
tanh x= sinh x 
cosh x 


Inverse hyperbolische Funktionen: 


sinh“'x= In| x-+ (x?+ 1)"/2] 


cosh"!x= In[x+(x2- 1)1/2] x>1 


antrat | | x2<] 


@ Das Programm des Moduls belegt das Register 00. 


® Das Drucker-Einschalt-Flag (Flag 21) wird vom Modul-Programm nicht 
gesetzt. 


Bemerkungen: 


Hyperbolische Funktionen 45 


Umfang: 001 


am Tan Trmim 


Für hyperbolische Funktionen weiter 
mit Schritt 2; für inverse hyperbolische 
Funktionen weiter mit Schritt 3. 


HYPERBOLISCHE FUNKTIONEN 


Argument eingeben und berechnen 

© Sinus hyperbolicus 

© Kosinus hyperbolicus 

© Tangens hyperbolicus 

INVERSE HYPERBOLISCHE FUNKTIONEN 


Argument eingeben und berechnen 
© inverser Sinus hyperbolicus 

@ inverser Kosinus hyperbol. 
© inverser Tangens hyperbol. 


Beispiel ! 
Berechne die folgenden hyperbolischen Funktionen: 

sinh2.5, cosh3.2, tanh1.9 
Tasten: Anzeige: 
SIZE 001 
2.5 SINH 6.0502 (sinh 2.5) 
3.2 COSH 12.2866 (cosh 3.2) 
1.9 TANH 0.9562 (tanh 1.9) 
Beispiel 2: 


Berechne die folgenden inversen hyperbolischen Funktionen: 


sinh’! 2.4, cosh’! 90, tanh”! -0.65 


Tasten: Anzeige: 

2.4 ASINH 1.6094 (sinh"! 2.4) 
90 ACOSH 5.1929 (cosh”! 90) 
.65[CHS) 


ATANH -0.7753 (tanh"' —0.65) 
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DREIECKSBESTIMMUNGEN 


Diese Programme können zur Bestimmung von Fläche, Seitenlänge (S,, S,, S;) 
und Winkel (A,, A», As) ebener Dreiecke eingesetzt werden. 


A, 


Es werden einfach die drei bekannten Werte eingegeben und das entsprechende 
Programm gewählt. Der Rechner gibt Fläche, Seitenlängen und Winkel aus. Die 
Reihenfolge der Ausgabe hängt von der Reihenfolge der Eingabe ab. Werden 
die Eingabewerte nach dem Uhrzeigersinn um das Dreieck gewählt, so werden die 
errechneten Werte auch im Uhrzeigersinn ausgegeben. Die Reiheniolge ist: 


erste eingegebene Seite (S}) 
anliegender Winkel (A,) 


anliegende Seite (S,) 
anliegender Winkel (A,) 
anliegende Seite (S;) 
anliegender Winkel (A;) 
Fläche 

Gleichungen: 


Alle Dreiecksseiten (S,, S,, S;) bekannt: 


P(P — S,) 


A = -1 
>= 2.008 5,5, 


wobei P= (Sı + S, + S3)/2 


E PP-S 
Aum2conı , ET 


A, = cos”! (-cos(A; + A3)) 


Dreiecksbestimmungen 


A;, Sı und A, bekannt: Zwei Winkel und die eingeschlossene Seite. 
A, = cos! (-cos(A, + A,)) 


=S, ni 
S; = S, cos A, + $S, cos Ay 

S], A, und A, bekannt: Eine Seite und die beiden folgenden Winkel. 
As = cos”! (-cos(A, + A3)) 


Damit ist das Problem auf die A;, S,, A, Lösung reduziert. 
S], A, und S, bekannt: Zwei Seiten und der eingeschlossene Winkel. 


Sz = \ S,2 + 8,” -2 SIEH cos A, 


Damit ist das Problem auf die S,, S,, S; Lösung reduziert, 
Sı,S, und A, bekannt: Zwei Seiten und der anliegende Winkel. 


. * 
A; = sin"! E sin | 
Sı 


A, = cos-![ -cos(A, Hr A;)] 


Damit ist das Problem auf die A;, S,, A, Lösung reduziert. 


Area = 1/2S, S, sin A, 
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*Beachte: Es existieren zwei Lösungen für den Fall, daß Sz größer als Sı und Az ungleich 90° ist. 


Beide möglichen Lösungen werden bestimmt. 


48 Dreiecksbestimmungen 


Anmerkungen: 
@ Das Programm belegt die Register Ryy-Ror. 
@ Winkelangaben müssen dem Winkelmodus des Rechners entsprechen. 


@ Es ist zu beachten, daß das Dreieck in der Programmbeschreibung nicht der 
Standardnotation für Dreiecke entspricht, z. B. liegt A, der Seite S, nicht gegen- 
über. 


@ Winkel müssen als Dezimalzahl angegeben werden. Die[HR ] Konversion kann 
dazu verwendet werden, Grad, Minuten und Sekunden in dezimale Winkel- 
angaben umzuformen. 


@® Die numerische Genauigkeit kann für sehr kleine Winkel ungenügend sein. 


entsprechenden Winkelmodus einstellen 


anwendbaren Fall der untenstehenden 
Liste wählen und die angezeigten Werte 
eingeben: 


alle Seiten bekannt 


Zwei Winkel und die eingeschlossene 
Seite bekannt 


Zwei Winkel und anliegende Seite 
bekannt 


Zwei Seiten und der eingeschlossene 
Winkel bekannt 


Zwei Seiten und der anliegende 
Winkel bekannt 
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Nach Schritt.2 können die Werte für 
Seiten und Winkel durch sukzessives 
Drücken von angezeigt werden. 
Die letzte Angabe ist die Dreiecksfläche. 
Für den letzten Fall (SSA) gibt es zwei 


mögliche Lösungen, beide werden 
ausgegeben 


Beispiel !: 
Gesucht sind die Winkel (in Grad) und die Fläche des folgenden Dreiecks: 


2.75 


Tasten: Anzeige: 


58 s1-3 

2 Ss2-? 

1 s3=? 
2.75[R2) S1=2.0000 
A1=129.8384 
un 
S3=2.7500 

A3=16,2136 

R/S 


AREA=0.7679 
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Beispiel 2: 

Ein Landvermesser muß die Fläche und die Ausdehnung einer dreieckigen 
Parzelle bestimmen. Vom Punkt A werden die Abstände zu B und C mit einem 
elektronischen Längenmeßgerät bestimmt. Die Winkel zwischen AB und AC 
werden ebenfalls gemessen. Gesucht sind Fläche und die restlichen Seitenlängen. 


Dies ist ein Seite-Winkel-Seite Problem mit: 


S =171.63, A,=98° 12° und S,=297.35 


Tastenfolge: Anzeige: 
SAS s1=? 
171.63 [R/S)  Al=? 
98.12 HR 

S2=? 

297.35 S1=171.6300 
R/S -A1=98.2000 
S2=297.3500 

A2=27.8270 
S3=363.9118 
R/S A3=53.9730 


R/S AREA=25,256.2094 
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Beispiel 3: 

Ein Pilot möchte genau nördlich fliegen. Die Windangabe ist 25 Knoten und 
077°. Da Windrichtungen diejenige Richtung sind, aus derder Wind kommt, wird 
der Windwinkel in 077° + 180° oder 257° umgerechnet. Die wahre Eigen- 
geschwindigkeit des Flugzeug ist 140 Knoten. Welcher Steuerkurs muß - zur 
Kompensation der Abdrift durch den Wind - geflogen werden? Wie hoch ist die 
Geschwindigkeit über Grund? 

Steuerkurs: HDG (HeaDinG) wahre Eigengeschwindigkeit: TAS (True Air 
Speed) 

Geschwindigkeit über Grund: GS (Ground Speed) Windgeschwindigkeit: WV 


N 


Zieht man die Windrichtung von 180° ab (man erhält dabei einen Winkel von 
103°), reduziert sich das Problem auf eine Sı, S2, A» Lösung. 


Tastenfolge: Anzeige: 
SSA s1=? 
140 [R/S) S2=? 
25 A2=? 
103 S1=140.0000 (TAS) 
A1=66.9798 
S2=25.0000 (Wind velocity) 
R/S A2=103.0000 
S3=132.2407 (GS) 
R/S A3=10.0202 (HDG) 
R/S AREA=1,610.6428 


Somit muß der Pilot einen Steuerkurs von 10.02° (NNO) fliegen. Seine Ge- 
schwindigkeit über Grund beträgt 132.24 Knoten. 


52 
KOORDINATENTRANSFORMATION 


xfyotzu atbtctO x*’ytz-P 


Dieses Programm bietet 2- und 3-dimensionale Koordinaten-Translation 
und/oder Rotation. 


Für den 2-dimensionalen Fall sind die Koordinaten des Ursprungs des achsenver- 
schobenen System und der Rotationswinkel relativ zum Ausgangssystem einzu- 
geben, danach sind die neuen Achsen zu spezifizieren. (Ursprungskoordinaten: 
(X0,Y0), Rotationswinkel: 9). Diese Größen werden mit der Taste [AJeingegeben. 
Danach können durch Drücken von [c] Punkte des ursprünglichen Systems mit 
(x,y) eingegeben und in das neue achsenverschobene und gedrehte System mit 
(x’,y’) abgebildet werden. Punkte des neuen (x’,y’)-Systems können durch Drük- 
ken von] in das ursprüngliche (x,y)-System abgebildet werden. 


YSTEM 
Rn - 
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Der 3-dimensionale Fall ist dem 2-dimensionalen analog. Der einzige wichtige 
Unterschied betrifft die Angaben zur Drehung. Die Rotationsachse geht durch 
den verschobenen Ursprung (X, Yo, Zu) und ist parallel zu einem beliebigen Orts- 
vektor (ai,bj, ck). Das Vorzeichen des Rotationswinkels (©) wird nach der Rechte- 
Hand-Regel und der Richtung des Rotationsvektors bestimmt. Beispielsweise in 
dem Spezialfall der 2-dimensionalen Rotation (Drehung nur in der (x,y)-Ebene) 
kann dies unter Verwendung eines Ortsvektors von (0,0, 1) und einem positiven 
Drehwinkel - für Drehung entgegen dem Uhrzeigersinn- bewerkstelligt werden. 
Die Koordinaten des verschobenen Ursprungs (X, Yo, Z0) werden mit m [AJein- 
gegeben. Der Ortsvektor und der Drehwinkel werden mit m eingegeben. 
Umwandlungen von ursprünglichen Systemen (x,y,z) in das neue System mit 
(x’,y’,z’)werden mit ı [c] initiiert, während der umgekehrte Vorgang mit ıı [E] 
erreicht wird. 


x fı m, Nı XT7% 
y u. f: m N y7”Yyo 
zZ Is Ma; N3 y Aa 7 \ 


' 
y = m; Ma Mg; y T Yo 
zZ nı Na Nn3 z’ Zy 
wobei: 


bı % ls | Ja®tı-cos#)+cos® ab(l-cos#)-csin ac(1-cos #)+bsin® 
m, m, ms[=|ba(l-cos#)+csin@ b?(1-cos #)+cos® be(l-cos#)-asin® 
N; ng nz | [call-cos#)-bsin® cb(1-cos#)*tasin 9 c*(1-cos@)+cos 9 


Zweidimensionale Transformation werden als Spezialfall der Dreidimensionalen 
Transformation behandelt; dabei wird (a,b,c) auf (0,0,1) gesetzt. 


Anmerkungen: 

@® Für reine Translation ist 9 = 0.0 zu setzen. 

® für reine Rotation sind X,, yo und Zu gleich Null zu setzen. 
© das Programm belegt die Register 00-25 


eg: 026 


| Nummer] Anweisung —___ | Engel rundum _ | Anzei 


Programm initialisieren und Maske 
auflegen. 

für den 2-dimensionalen Fall: 
weiter mit Schritt 3 

für den 3-dimensionalen Fall: 
weiter mit Schritt 6 


TRANS | 0.0000 
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Ursprung und Rotationswinkel des neuen 
Koordinatensystems eingeben 


sul 


Transformieren der Koordinaten 
des alten Systems in das neue 


"fl DB] 


oder 


aus dem neuen System in das 
alte 


7 


für ein weiteres Koordinatenpaar weiter 
mit Schritt 4. Bei einem neuen 
Koordinatensystem weiter mit Schritt 3. 


Ursprung des verschobenen Koordinaten- 
systems eingeben 


= 


und 


Rotationsvektor und -Winkel 
eingeben 


5 


Va+r+02 


2 
®) 


Koordinten des ursprünglichen Systems 
in das neue System umwandeln, 


iR 


22 
— > 
33 


oder 


aus dem translozierten und rotierten 
System in das ursprüngliche 


Zeit]: 


für ein neues Koordinatenpaar weiter 
mit Schritt 7. 

für eine neue 3-dimensionale Transfor- 
mation: weiter mit Schritt & (sowohl 
(Xo, Yo: Zu) als auch (a,b,c,@) können 
unabhängig voneinander geändert 
werden. 
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Beispiel 1: 


Gegeben seien die beiden untenstehenden K.ooordinatensysteme (x,y) und (x’,y”): 


U} 
© Pı(-9, 7) 


P:(6, 8) 


\ 
u = 27° 
\ 


P& (2.7, -3.6) 
1} 


Die Punkte P,, P, und P; sollen in die äquivalenten Koordinaten des (x’,y”)- 
Systems überführt werden. Außerdem soll P, in das (x,y)-System überführt 


werden. (Es werden Winkelgrade verwendet.) 


Tastenfolge: Anzeige: 


SIZE [ALPHA] 025 
TRANS 0.0000 
7 a[cHs) 


27[a) 1.0000 
9[cHS) 7) -9.2622 
17.0649 
5[CHS) a[cHs)[C) -10.6921 
5.4479 
6 [Erere) sc) 4.5569 

R/S 11.1461 
2.7 3.6[cHs)[E) 11.0401 

R/S -5.9818 
Beispiel 2: 


(X) 
(y 1) 
(X,') 
(y2) 
(X) 
(y3) 
(X4) 
(y) 


Ein 3-dimensionales K.oordinatensystem wird in den Punkt (2.45, 4.00,4.25) ver- 
schoben. Danach wird eine Drehung um 62.5° um die Achse (0,-1,-1) ausgeführt. 
Im ursprünglichen System hat ein bestimmter Punkt die Koordinaten (3.9, 2.1, 
7.0). Was sind die Koordinaten dieses Punktes im neuen System? 


56 Koordinatentransformation 


Tastenfolge: Anzeige: 

TRANS 

2.45 4 (EnteRe) 

4.25 (A) 2.4500 

0 1 

1 62.5 1.4142 

3.9 (TER?) 2.1 

7 3.5861 (x) 
0.2609 (y’) 
0.5891 (z') 


In dem neuen (verschobenen und gedrehten) Achsensystem hat der Punkt die 
Koordinaten (1, 1,1). Was waren seine Koordinaten im ursprünglichen System 


Tastenfolge: Anzeige: 

1 1 

ı@8 (E) 2.9117 (x) 
4.3728 (y) 


5.8772 (z) 
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ANHANG A 
PROGRAMM-DATEN 


1 
17 
20-00 &7-00 05 LONEUNOJSUENUSJEUIPIOOY 
IRA 10-00 gr uaßunuysasaqsyaaıaıq 
00 il uauonyung ayasıjoqsadAH 
930/QVY IZ +0-00 65 uauoneladg axajdıoy 
67 
14 
vXl4 °0-10 97-00 08 askjeuy-IalınoJ 
[R4 
0 10-00 cE usßunyslsjäjenussayig 
13 
LZ 10-00 IR uonesbalu] aydsuauunn 
1% 
&0-20 u3WouAjod UOA ayaM 
YXl4 00 °7-00 v8 -suoyun, pun uajja}sjjnN 
07-17 
00 90-00 137 (X)} any uaßunso] 
67 
GL 
37-17 (uejdsßunßajaq 
01-20 -J2]sıbay 'IBA) BEI eigabjexıneyy 
SsnaoıN LvVNu0oJ u31S1934 AdOD UNI 
-134NIM -J91IZNV savVu4 N3lva IHVZU3LSIDO34Y NNVHDOHd 


EEE EEE EEE SEIEN EEEEENEEENEEENNSEENEEEEREEREEREEE EEE nn nn 


"UBPJIBOM JYS18119 PL +l+(}-I)N = 3553Jpe13}s1B34 :J9WJOJ 13p BunpusMmuaA Jajun uuey 


fi 


"e JUSWE|gZ SEpaf "Uaysıadsaß asıamugjlaz Is} XuJeW Old 


882-533 vzz-Liz 0OL2-S1 8E2-00 683 E 125 

602-261 961-F8l E8L-Sl 6023-00 0L2 E el 

OLr-Eel Sal-Lzl 0L1-6S1 851.-S1 28,-00 esl [4 [43 

9pL-E1 ZSL-Zpl 9P1-9E1 Ser-St ZSL-00 851 [4 rt 

vzr-el vel-Ssaı var-Ssit yLı-Ssı vEL-00 gel ° OL 
volr-ei ELL-SOL v01,-96 g6 St ELL-00 123° t 6 
98 -£L v6 -28 98 -62 82 -Sı v6 -00 S6 [ 8 
02 Er ZZ -VZ 02 -#9 e9 SI 22 -00 8 L L 
99 -£1 9 -1S 99 -15 05 -Sı <9 -00 &9 0 9 
vv el 6r -Sy vr -Or 6E -Sı 6r -00 05 0 S 
ve el BE -SE ve -IE 0& -Sı ge -00 6E 0 r 
93 Ei 62 "25 9 v3 e2 -Sı 62 -00 0E 0 e 
08 -£i 2 17 03 -61 8 -Sı cz -00 ez 0 [4 
9, -EL 2b = 9, -9L Sr -Sı ZL -00 8 (6) L 

YL+N+zN YL +NZ+zN YL+N+zN YL+zN FL +NZ+zN 
or eL 0} SL+N+ZN 0} SL+zN 0, SL 0,00 SL+NZ+ZN (tr) INI N 
(w) 
(wvy) aınpow-sBunu 
uayeds ı3p S}OAId J3P xuyeW J3p je10} (4) 191516534 -ayamIslaysıads 
494s1534 9SSFIPEPUZ pun -sbugjuy Iyezuy Iyezuy uolsuawig 


N3ZIULVIN UNI NVIIdSONNDFT3HUFIHD134S 
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ANHANG B 
UNTERPROGRANIME 
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HP-41C 
Mathematik-Paket | 
KURZBESCHREIBUNGEN 
DER PROGRAMME 


MATRIXALGEBRA 
Programm initialisieren: xEQ] MATRIX 


Anzahl (N) von Zeilen/Spalten eingeben, (N=14); 
drücken. 


Speicherumfang eingeben und weiter. 


Matrixelemente Ajj zeilenweise eingeben, [R/S 
drücken. 


Diesen Schritt wiederholen, bis alle Elemente 
eingegeben sind. 


Matrix mit [xeo] VIEW ausgeben. 
Matrix mit [xeo] EDIT ändern: 


Zeile i und Spalte j des neuen Elementes mit i 


ENTER4] j eingeben, drücken. 
Neuen Wert Aij eingeben, drücken. 
Änderungsroutine EDIT mit [Rs] [Rrs] beenden. 


Printed in SGP c Hewlett-Packard 
1979 


Determinante berechnen: [xeo] DEI. 

Inverse bestimmen: [xeo] INV. 

Ergebnisse spaltenweise ausgeben: [RS 
Lösen simultaner Gleichungen: SIMEQ. 


Spaltenvektor B eingeben: Bı 
restliche Werte von BB ein- und Ergebnis aus- 


geben: [RıS. 
Spaltenvektor ausgeben: [xcEo| VCOL. 


LÖSUNGEN VON fx) =O IN EINEM INTERVALL 
Mindestumfang: 007 

Funktion unter gewünschtem/GTo]| [e] [®] 

Label eingeben. 


RTN 
Programm initialisieren: xeo] SOLVE 
Funktionsnamen eingeben: [Rs | drücken. 
Sollen 2 Anfangsschätzungen eingegeben 
werden, deren erste eingeben: 


andernfalls leere Eingabe mit [rvS. 
Zweite Anfangsschätzung mit eingeben. 


NULLSTELLEN UND FUNKTIONSWERTE 
VON POLYNOMEN 


MINDESTUMFANG: 023 


Programm initialisieren: xEo] POLY 
Polynomgrad n (n=2,3,4,5) eingeben: 


Koeffizienten aj eingeben: 
Schritt wiederholen bis ROOTS? in der Anzeige. 


Um Lösungen zu bestimmen: drücken. 


Um Funktionswert zu berechnen, mit nein (N) 
antworten, drücken. 


x eingeben, drücken um f(x) auszugeben. 


Für ein neues Polynom gle'chen Grades, die | 
Koeffizienten in Roo-Ro4 ändern und [xcEo 
ROOTS. 


NUMERISCHE INTEGRATION 
Mindestumfang: 008 


Diskreter Fall | 
Programm initialisieren: Xeo] INTG. 


Abszissenintervallgröße eingeben: [a]drücken. 


Funktionswert für xj eingeben; [B] drücken. 
Für j=0,1,..., n wiederholen. 


Fläche nach der Trapezregel: [C], nach 
SIMPSON berechnen: [D|. 


Explizite Funktionen 


Funktion unter gewünschter [e] [e] 
Marke eingeben. 


LB 


T 


I 
2 rm 


Programm initialisieren: [xEo| INTG. 


Intervallgrenzen a und b eingeben: 
a [ENTER bE A) 


Anzahl der Teilintervalle eingeben und Fläche 
nach SIMPSON berechnen: 9 [B]. 


Funktionsnamen eingeben, drücken. 


Sollen a,b, oder n geändert werden, so ist mit 
dem entsprechenden Schritt fortzufahren. 


DIFFERENTIALGLEICHUNGEN 
Mindestumfang: 008 
Funktion eingeben: [(GTO] 


T 


2 m 0) 
ERS 
Q 
Q 


Programm initialisieren: [xEo| DIFEQ 
Funktionsnamen eingeben: drücken. 


Ordnung der Differentialgleichung (1 oder 2) 
eingeben, drücken. 


Schrittweite eingeben: 
Anfangswert für x eingeben: 
Anfangswert für y eingeben: 


Für Gleichungen 2. Ordnung Anfangswert 
für y’ eingeben: [RıS). 
Aufeinanderfolgende Werte fürx und y 
durch Drücken von ausgeben. 


FOURIER-ANALYSE 

Mindestumfang: 008 

Programm initialisieren: XEo] FOUR 

Anzahl Stichproben in einer Periode eingeben; 
drücken. 

Anzahl der gewünschten Frequenzen eingeben; 
drücken. 

Yn mit n=1,...,N eingeben; [rs] drücken. 
Schritt wiederholen bis RECT? angezeigt wird. 
Zur Ausgabe in kartesischen Koordinaten: [rıs]. 
Ausgabe der Koeffizienten in Polarform, 

N eingeben; drücken. 


Aufeinanderfolgende Koeffizienten mit 
ausgeben. 


Um den Wert von f (t) zu berechnen, USER-Modus 
einschalten, t eingeben und [E] drücken. 


KOMPLEXE 
OPERATIONEN 
Mindestumfang: 005 


ARITHMETIK 
Erste komplexe Zahl (x1 + iyı) eingeben: 


y1 x1 

Zweite komplexe Zahl (x2 + iy2) eingeben: 
y2 [ENTER+] .X2 

+ C+ 

eg C- 

x KEOol Cx 

= ea] C/ 


FUNKTIONEN 
z y1 x4 MAGZ 
2° y1 x1 CINV [RiS] 


z"  yı [ENTERS] X1 [ENTERA] n 
KXEO] ZTN [Rs] 

zn  yı [ENTERA] xı [ENTERS] n 
XxeEoO] Z T1/N [rs] 

e y1 [ENTER$) %1 XEo] eTZ [rs] 


Inz yi [ENTERt] X1 xEo] LNZ [rıs] 


Eu y1 [ENTERt x1 [ENTERt a 
xEQl aTZ [RıS] 


al092 y4 [ENTERR] x1 [ENTERR) a 


XEO] LOGZ 
HEWLETT ä 
@ Be 00041-90066 


2. 3 x2 y1 
x1 XEQ] zZ /W [Rs] 
zIw KEOl z 11IW RS 


sinz Yy1 |ENTER#| x [xEQo| SINZ [RıS 
coSZ Yy1 [ENTERA4| %X1 (XEQ] COSZ [rıs 
tanz y1 [ENTER4| X1 (xEo| TANZ [RS | 


HYPERBOLISCHE FUNKTIONEN 
Mindestumfang: 001 


sinn x x[xeEo| SINH 
cosh x x[xcEo] COSH 
tanh x x [xeEol TANH 
sinn’! x x XEo] ASINH 
cosh !x x KEG] ACOSH 
tanh’!x x eo] ATANH 


DREIECKSBERECHNUNGEN 
Mindestumfang: 008 | 
Alle Seitenlängen bekannt: [xEo] SSS 


Zwei Winkel und die eingeschlossene Seite 
bekannt: [xEQ| ASA 


Zwei Winkel und die anliegende Seite bekannt: 
XEQ| SAA 

Zwei Seiten und der eingeschlossene Winkel 
bekannt: [xEQ| SAS 


Zwei Seiten und der anliegende Winkel bekannt: 
XEQ| SSA 


Sind alle Fragen beantwortet, werden die 
Ergebnisse mit aufeinanderfolgenden aus- 
gegeben. 


KOORDINATENTRANSFORMATION 
Mindestumfang: 026 
Programm initialisieren: XEQ| TRANS 


2-DIMENSIONAL 
Ursprungskoordinaten und Drehwinkel eingeben: 


Xo [ENTER] Yo [ENTERt| 9A} 


Koordinaten in verschobene und gedrehte über- 
führen: x |ENTER#| y [c] [rs]. 


Koordinaten rücktransformieren: 


x |ÄENTER#| y’ [E] [rıs]. 


3-DIMENSIONAL 


Ursprungskoordinaten eingeben: 
Xo [ENTER Yo [ENTERt] Zo E [A| 
Drehvektor und -winkel eingeben: 


a ENTER#| b |ENTER#| c |ENTERt| 9 W lB] 


Koordinaten ins verschobene und rotierte System 
transformieren: 


x [ENTER#] y [ENTER#] z m [C] [Rs] [Rs 


Koordinaten rücktransformieren: 


x [ENTERt] y’ [ENTER ZU ERS] 


HP-41C 


Math Pac 
Addendum Card 


In the program “TRANS”, clear Flag 01 before 
attempting to transform a point from the 
translated system back to the original system. If 
Flag 01 is not cleared, the results will be in error. 
This applies to both 2 and 3 dimensional cases. 


The program “POLY” may calculate incorrect 
roots for fourth and fifth degree equations. To 
determine whether a calculated root is correct, 
rerun the Poly program and follow the in- 
structions on page 22 for evaluating a function. 
When the program asks “X = ?”, input the root in 
question. Iftheresult is a very small number close 
to zero, then the root is correct. For fifth degree 
polynomials, only the second real root need be 
evaluated. For fourth degree polynomials, only 
the first real root need be evaluated. 
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